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Abstract 

A detailed NNLO QCD analysis of new CCFR data on the xF^ structure func- 
tion including the target mass, higher twist and nuclear corrections was performed 
and parametrizations of the perturbative and power terms of the structure func- 
tion were constructed. The results of QCD analysis of the structure function were 
used to study the Q^— dependence of the Gross-Llewellyn Smith sum rule. The 
a^/vr— expansion oi ScLsiQ'^) was studied and parameters of the expansion were found 
to be si = 2.74 ± 0.01, S2 = -2.22 ± 0.23, S3 = -7.86 ± 1.74 in a good agreement 
with the perturbative QCD predictions for the Gross-Llewellyn Smith sum rule in the 
next-to-leading and next-to-next-to- leading order. 
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1 Introduction 

The progress of perturbative Quantum ChromoDynamics (QCD) in the description of 
the high energy physics of strong interactions is considerable Recent experimental data 
on the structure function of neutrino deep-inelastic scattering obtained at the Fermilab 
Tevatron p| provide a good possibility to precisely verify the QCD predictions for scaling 
violation with experiment data. The QCD predictions for the structure functions evolution 
of deep-inelastic scattering (DIS) are calculated now up to the next-to-next-to-leading 
order (NNLO) of a perturbative theory |P, ||, |^. The method of comparison of 3-loop 
QCD predictions with SF experimental data has been developed in [0, ||, ^ based on the 



Jacobi polynomial structure function expansion [|10[. It is well known that beyond the 
perturbative QCD there are other effects (higher twist effects, nuclear corrections, target 
mass corrections etc.) to be included into the joint QCD analysis of SF. 

In the paper, we present the results of NNLO QCD analysis of the data on the xF^ 
structure function obtained by the CCFR Collaboration with taking in account of the target 
mass, higher twist and nuclear corrections. The Q^-evolution of SF is studied and the 
parametrizations of perturbative (leading twist) and power terms of the structure function 
are constructed. The results of our NNLO QCD analysis of the structure function are in 
good qualitative agreement with NNLO theoretical predictions for the Q^-evolution of the 
Gross-Llewellyn Smith sum rule ^ S'^l°^{Q'^) = 3 ■ [1 - a^/vr - 3.25 ■ (as/vr)^] 11. 

In Section 2, the method of NNLO QCD analysis of SF based on the SF Jacobi polyno- 
mial expansion including the target mass, higher twist and nuclear corrections is described. 
The results of NNLO QCD analysis of SF are presented in Section 3. In Section 4, the as/vr- 
expansion of the Gross-Llewellyn Smith sum rule is considered and expansion parameters 
are found. 



2 Method of QCD analysis 

2.1 Jacobi Polynomial Expansion Method 

We use, for the QCD analysis, the Jacobi polynomial expansion method proposed in [PH]]. 
It was developed in |]13[ and applied for the 3-loop order of perturbative QCD (pQCD) to 
fit F2 and xFg data g |, |l5l. 

Following the method [1^ , we can write the structure function xF^ in the form: 

^ max 

xFf^^ix, Q') = x-il - xf Y: 0n'"(^) E /3)M|'3CD ^ 2, Q') , (1) 

n=0 j=0 

where 9"^(x) is a set of Jacobi polynomials and c^{a,P) are coefficients of the series of 
G^'^(x) in powers of x : 

n 

e:'^{x)=j:cf{a,(3)x^. (2) 

j=0 

The - evolution of the moments M^^^^ [N, Q^) is given by the well-known perturba- 
tive QCD HI, |13 formula: 



2 



rPQCD 



Mr {N,Q 



«5 m 



as 



H^{Ql,Q')Mf^''{N,Ql), AT = 2,3,... (3) 



Here s(Q^) is the next-to-next-to-leading order strong interaction constant, 7^''^"^ are 
the nonsinglet leading order anomalous dimensions, and the factor Hjy (Qq, Q^) contains 
next- and next-to-next-to-leading order QCD corrections to the coefficient functions and 
anomalous dimensions i and is constructed in accordance with 0, ^ based on theoretical 
results of l-g. 

The unknown coefficients M^{N, Qq) in (|^) could be parametrised as Mellin moments of 
some function: 

Mf^^{N,Ql) = (ixx^-^aix^^il-xY^il + a^x), N = 2,3,... (4) 
Jo 

Here coefficients Oj should be found by the QCD fit of experimental data. 

The target mass corrections (TMC) are included into our ffis through the Nachtmann 
moments [jl9| of the SFs: 

MPQC^{Q^) = ['dxC^'/{x^)F,{x,Q'){l + {n + l)V)/{n + 2), (5) 
Jo 



where ^ = 2x/ (1 + V), V = + AM^^^ix'^ / Q'^ and Mnud is the mass of a nucleon. We are 
taking into account the order 0{M^^^i/Q^) corrections: 

MniQ') = M„^«^^(g2) + ^^I^^MSf^(Q^) (6) 

where are the Mellin moments of measured xF^ SF. 

2.2 Higher Twist Correction 



We take into account the HT contribution following the method of and |T^, [T^, y]. 
To extract the HT contribution, the SF is parametrized as follows: 

xF3(x,g2)=xFf«^^(x,Q2) + ^, (7) 

where the dependence of the ffist term in the r.h.s is determined by perturbative QCD. 

The function h{x) as well as the parameters ai — 04 and scale parameter A should be 
determined by fitting the experimental data. 

2.3 Nuclear Correction 



We have used the covariant approach in the light-cone variables PT| , P2[ to estimate the 
ratio of structure functions 

j^D/N _ F^{x,Q'^) 



^For reviews and references on higher order QCD resuhs, seepS] 
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and to perform the joint NNLO QCD analysis of the data on the structure function 



(for the NLO analysis, see p3 



We would like to remind that the ratio Rp^^ describes the influence of nuclear medium on 
the structure of a free nucleon in the process. We use the approximation that Rp^^ = Rp^^^ . 
It gives lower estimation on the effect of nucleon Fermi motion in a heavy nucleus. 

The covariant approach in the light-cone variables is based on the relativistic deuteron 
wave function (RDWF) with one nucleon on mass shell. The RDWF depending on one 
variable, the virtuality of nucleon k'^{x, k±), can be expressed via the Dpn vertex function 
Va{x^ k±). The deep-inelastic structure function of neutrino-deuteron scattering in the 
approach can be written as follows 

F3^(a, Q2) = dx d^k^ A{x, k^) ■ F.^{a/x, Q^). (9) 

J a 

The nucleon SF is defined as Fj^ = {F^^ + if^)/2, a = -q^/2{pq). The function 
A(a;, kj_) describes the left (right)-helicity distribution for an active nucleon (antinucleon) 
that carries away the fraction of deuteron momentum x = ki^/p^ and transverse momentum 
k±. It is expressed via the RDWF V'a(^i) as follows 

A(x, k^) oc Spirih) -{m + k)- ^^(A;i) ■ q ■ a^'^ ■ ■ p^^ " W?V}, (10) 

is) 

where p^,^ is the polarization density matrix for an unpolarized deuteron. Note that in 
the approach the use of the distribution function A(x, ki_) includes not only usual S- and 
D-wave components of the deuteron but also a P-wave component. The latter describes 
the contribution of A^iV-pair production. The contribution of this mechanism is small over 
a momentum range {x < 1). 

The nuclear effect in a deuteron for the u + D ^ fi~ + X process has been estimated in 
p3| . It has been found that the ratio Rp^^ is practically independent of the parametrization 
of parton distributions |25|, |26| and the nucleon SF over a wide kinematic range of 



X = 10 ^ — 0.7, = 1 — 500 {GeV/c)"^. The curve has an oscillatory feature and cross-over 



point Xq : Rp^'^ {xo) = 1, Xq — 0.03. 



Thus, the obtained results give evidence that the function Rp^^ is defined by the struc- 
ture of the relativistic deuteron wave function and can be used to extract the nucleon SF 
F^ from the experimentally known deuteron one 

F3^(x, Q') = ■ Fi^ix, Q'). (11) 

The performed analysis of the nuclear correction for the nucleon SF also allows one to 



consider the influence of the nuclear effect on the GLS sum rule 11 



Sgls= C F^{x)dx. (12) 
Jo 

We have used the result on the R^^^ ratio to study the x and Q^-dependences of the 
GLS integral 

Sgls{x,Q^)= f F^{y,Q^)dy. (13) 

J X 

It has been shown in ||23| that the nuclear effect of Fermi motion is very important for 



the NLO QCD analysis of the structure function F^ and verification of the Gross-Llewellyn 
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Smith sum rule over a wide region of = 3 — 500 [GeVjc)^. Therefore in the paper, we 
take into account the nuclear effect within the joint NNLO QCD analysis of the structure 
function and of the Gross-Llewellyn Smith sum rule. 



3 NNLO QCD analysis of Structure function xF^ 

In this section, we perform the QCD analysis of the xF^ experimental data taking 
into account the 3- loop QCD, higher twist and nuclear corrections. We consider as a first 
approximation that i??/^ = RT = R- 

xFi'ix, Q') = ai(g2)a;'^2(Q^)(i _ xr^Q'\l + a,{Q')x) + (14) 

Constants h{xi) (one per x-bin) parametrize the HT x-dependence. The points Xi are chosen 
in accordance with 0: 0.0075 - 0.75. 

The values of the parameters ai{Ql) (i = 1 — 4) , scale parameter A and constants h{xi) 
have been determined by fitting the whole set of xF-^ data (116 experimental points) 
for different values of Ql in the kinematic region: 2 GeV^ < < 200 GeV"^. □ Only 
statistical errors are taken into account. The parameter is found to be about 105 for 116 
experimental points. 

Figure l(a-d) shows the parameters SF Oi at the points = 1.5 - 200 (GeV/c)'^. The 
Ajjg parameter was found in the interval (210 — 250) MeV with statistical error about 
±20 MeV. 

Figure 2 demonstrates the x-dependence of the HT contribution h{xi). 
We use the result of our NNLO QCD analysis presented in figs. 1 and 2 to obtain the 
parametrisation of xF^ SF in form (15). 

The Q^-dependence of the coefficients is parametrized as follows: 

a,{Q') = j2crz^, z = \og{Q') (15) 

j=0 

and is shown in Fig.l by solid lines. The coefficients q are presented in Table 1. 
Table 1. Coefficients q for the parametrization of the function xF3{x,Q'^) 





ai 


^2 


0-3 


04 


Co 


3.1430 


6.6704e - 


1 


3.2753e + 


2.92880e + 


Cl 


1.6794 


4.0720e - 


3 


8.2294e - 1 


-1.65755e + 


C2 


-0.3896 


-7.7966e - 


3 


-1.3154e- 1 


3.09344e - 1 



Figure 1 shows the dependence of coefficients ai on Q^. The black circles are experimental 
points, lines are results of the fit. 

"'The number of flavors is taken to be equal to 4. 
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The parameter 02 slowly depends on Q"^ and corresponds to the theoretical estimation 

^theor ^ Q Qg 

Our results for 03 (Q^) could be compared to the well known Buras-Gaemers parametriza- 
tion ET 



a3 = Vi+m- 



(16) 



At the point Ql = 1.8 GeV^ it gives r/i = 2.6 and r]2 = 0.8. Our result is r]i = 3.48±0.03 
and rj2 = 1.18 ± 0.09. It is not far from the pQCD theoretical estimation based on the 
relation 

^a3(Q^) = A„,(g.) + 0(a|(Q^)) (17) 

and gives r^*'"^°'" = 0.64. 

The term describing the higher twist correction is written in the form 



i=0 

and is shown in Fig. 2 . The values of coefficients di are presented in Table 2. 
Table 2. Coefficients di for the parametrization of the function h{x) 



do 


di 


d2 


ds 


0.2329 


0.8060 


0.7308 


0.1842 



:i8) 



Thus, the NNLO QCD analysis of experimental data of the structure function 0:^3 has 
been performed and the parametrizations of the NNLO perturbative QCD, TMC, nuclear 
effect and higher twist corrections have been obtained. 



4 Experimental constraints on coefficients of a^-expansion 
of Gross-Llewellyn Smith sum rule 

In this section, we would like to show the status of the NNLO QCD analysis of the 
experimental data 0. 

The Q^-dependence of the parameters ai{Q^) allows us to study the behavior of Sgls{Q^) 
in a wide kinematic region of the momentum transfer squared. 

Figure 3 shows the dependence of Sgls{Q^) on as/'n' in thw 3-loop approximation. The 
black and open circles correspond to the NNLO QCD analysis with and without nuclear 
corrections, respectively. The statistical errors are about ±0.4 and are not presented at the 
figure. 

Figure 3 demonstrates the increase of Sgls{Q^) with decreasing as/ii- The result is in 
qualitative agreement with the Q^-dependence of the sum rule found in I^T] in the NLO 
QCD analysis without target mass corrections, higher twist and nuclear effect taking. A 
similar tendency in the NLO QCD approximation with target mass corrections and nuclear 



effect was found in For the estimation of the order 0{aj) and power corrections, see 



Figure 3 shows a considerable sensitivity of Sgls{Q^) to the nuclear correction. 
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The pQCD predictions for a^/vr-expansion of the Gross-Llewellyn Smith sum rule up to 
(as/vr)^ could be presented in the form 

SclsiQ') = Co ■ (1 + ci ■ {asM + C2 ■ {as/rrf). (19) 

The coefficients ci and C2 have been calculated in |jl2|: ci = —1, C2 = —3.25. 

The a5-dependence of Sqls presented in Fig. 3 could be parametrized by the parabola: 

SclsiQ') = So + Si- {asM + S2 ■ {as/nf. (20) 

This expansion allows us to directly compare the results of experimental data analysis with 
NNLO QCD calculations. The values of Sq, Si, S2 for the interval cts/vr < 0.09 are presented 
in Table 3. 

Table 3. The coefficients of the Gross-Llewellyn Smith integral Sgls{Q'^) expansion in 

as /it. 





NNLO+HT 


NNLO±HT±NC 


Theor. Q 


So 


2.86 ±0.01 


2.74 ±0.01 


3.00 


Si 


-4.93 ±0.32 


-2.22 ±0.23 


-3.00 


S2 


0.98 ±2.54 


-7.86 ± 1.74 


-9.75 



From Table 3 one can see a high sensitivity of parameters to the nuclear correction. The 
obtained results for the coefficients Si show that the consideration of the nuclear correction 
gives good agreement with the theoretical calculation of next-to-leading and next-to-next- 
leading order QCD corrections. 

5 Conclusion 

The detailed NNLO QCD analysis of the structure function xF^ of new CCFR data in- 
cluding the target mass, higher twist and nuclear corrections was performed. The parametriza- 
tions of perturbative and power terms of the structure function were constructed. The re- 
sults of QCD analysis of the structure function were used to study the Q^- dependence of 
the Gross-Llewellyn Smith sum rule. The as/vr-expansion of ScLsiQ"^) was exemined and 
expansion parameters si, S2, S3 were found. We would like to emphasize that the considera- 
tion of the nuclear correction allows us to achieve a good qualitative agreement between the 
results obtained from the experimental data by the NNLO QCD analysis and perturbative 
QCD predictions for the Gross-Llewellyn Smith sum rule in next-to-leading and next-to- 
next-to-leading order. 
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FIGURE CAPTIONS 



Figure 1. Dependence of the coefficients 01-04 on . The black circles are experimental 
points, lines are results of the fit. 

Figure 2. The function h{x) vs log{x) describing a higher twist correction. The black 
circles are experimental points, the line is the result of the fit. 

Figure 3. The cts/TT-expansion of the Gross-Llewellyn Smith sum rule Sgls{Q^)- Black 
and open circles are results of NNLO QCD analysis obtained with and without nuclear 
correction, respectively. 
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